NMR quantum simulation of localization effects induced by decoherence 
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The loss of coherence in quantum mechanical superposition states limits the time for which quan- 
tum information remains useful. Similarly, it limits the distance over which quantum information 
can be transmitted, resembling Anderson localization, where disorder causes quantum mechanical 
states to become localized. Here, we investigate in a nuclear spin-based quantum simulator, the 
localization of the size of spin clusters that are generated by a Hamiltonian driving the transmission 
of information, while a variable-strength perturbation counteracts the spreading. We find that the 
system reaches a dynamic equilibrium size, which decreases with the square of the perturbation 
strength. 
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PACS numbers: 03.65.Yz, 03.67.Ac, 72.15.Rn, 76.60.-k 



Introduction. — Quantum information processing has 
the potential of solving computational problems for 
which no efficient solution exists on classical computers 
[IHij]. Transfer and exchange of quantum information 
and quantum entanglement can be used for secure trans- 
mission of information 0-Q| • Realization of this potential 
for practical applications requires precise control of large 
quantum registers. However as the number of qubits in- 
creases, the quantum mechanical superposition states of 
the system become more fragile [||. This degradation of 
quantum superpositions, called decoherence [9j, is due to 
extra degrees of freedom (the environment) that interact 
with the system, and to imperfections of the gate opera- 
tions. Overcoming decoherence is clearly one of the key 
factors for implementing large scale quantum comput- 
ers. Several techniques have been proposed for this pur- 
pose, including dynamical decoupling [lOj, decoherence- 
free subspaces 11 1 , and quantum error correction fl2l . [l3j ]. 
These proposals have been tested on small systems of nu- 
clear spins , or trapped ions [ToT ] or spin model quan- 
tum memories fl6j |. 

Tests on larger systems, comprising hundreds or thou- 
sands of qubits, are more difficult. So far, the only phys- 
ical system that offered this possibility is nuclear mag- 
netic resonance (NMR) of dipolar coupled spins [H [l7| ■ 
Processes that transfer quantum information over large 
distances can also be studied in spin chains [5]. An exam- 
ple of such a linear spin system was studied by solid-state 
NMR [l8[ . These model systems do not allow addressing 
of individual qubits, but they allow one to study some 
aspects of decoherence and information transfer. In par- 
ticular, they can be used for studying the effect of the 
finite precision of experimental quantum gate operations 
on the transfer of quantum states: it was predicted that 
quantum information cannot be transmitted over arbi- 
trary distances, but that it will become localized [l9l - l23| . 

In this paper, we present the first experimental study 
trying to answer the following question: How far can 
quantum information be transmitted with quantum gate 
operations of finite precision? For this purpose, we use 
an NMR quantum simulator. Starting from individual, 



uncorrelated spins, we measure the build-up of clusters of 
correlated spins of increasing size. Introducing a pertur- 
bation to the Hamiltonian that generates these clusters, 
we find, that the size of the clusters reaches an upper 
bound. This upper bound appears to be a dynamic equi- 
librium: if the cluster size is initially larger than this 
equilibrium value, it decreases under the effect of the 
perturbed Hamiltonian, while the unperturbed Hamilto- 
nian leads to an increase. The equilibrium size decreases 
with increasing strength of the perturbation. 

Growth of spin clusters. — All the spins of the system 
are equivalent and they are in a strong magnetic field. In 
its Zeeman rotating- frame, the Hamiltonian of the spin 
system used for the quantum simulations is the high-field 
homonuclear dipolar interaction [24| 
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coupling constants. The quantum simulations start from 
the high-temperature thermal equilibrium [24] , p cx I z = 
Y^, P z - In this state, the spins are uncorrelated. 

We generate states with correlated spin clus- 
ters whose density operator terms are of the form 
7*. ..7^1^ (u,v,w — x,y,z), by letting the system evolve 
under the effective Hamiltonian 
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This Hamiltonian isprepared by means of a standard 
NMR sequence 2|| 26 1 shown in the upper part of Fig. 
[1] This Hamiltonian flips simultaneously two spins with 
the same orientation. Accordingly, the z-component of 
the magnetization M z changes by M — AM Z — ±2. At 
the same time, the number K of correlated spins changes 
by AK = ±1. 

To determine the average number of correlated spins, 
we use standard NMR techniques developed by Baum et 
al. . The technique relies on the fact that in a sys- 
tem of K spins, the number of transitions with a given M 
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Figure 1: (Color online) NMR sequence for the quantum sim- 
ulations. An unperturbed evolution is achieved when te = 0. 
The effective Hamiltonian Ho is generated by the sequence of 
7r/2 pulses shown in the upper part of the figure. 
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Figure 2: (Color online) Time evolution of the cluster size of 
correlated spins with the unperturbed Hamiltonian Ho- Dis- 
tributions of the squared amplitudes Am of density operator 
components as a function of the coherence order M are shown 
for two different cluster sizes. The latter are obtained from 
the half- width 2y / \n(2K) of the distribution function Am- 

shows a binomial distribution. For K 3> 1, the binomial 
distribution can be well approximated with a Gaussian of 
width oc \[K. To determine the effective size of the spin 
clusters in a given state, we decompose its density oper- 
ator p into components of coherence order M. They can 
be distinguished experimentally by rotating the system 
around the 2— axis: a rotation <f) z = e — by <j> changes 
the density operator to 

p{<t>)=hpfc 1 =Y.PMZ im , (3) 

M 

where pm contains all the elements of the density oper- 
ator involving coherences of order M. The terms with 
M = are zero quantum coherences and populations. 

If the system evolves under the Hamiltonian (|5J|, the 
cluster size increases indefinitely, as shown in Figure [5] 
The figure also shows two examples of pm distributions. 

This evolution can be reversed completely by changing 
the Hamiltonian from Ho to —"Ho- Experimentally, this 
is achieved by shifting the phase of all RF pulses by ±7r/2 
[25]. This indefinite growth of the cluster size, as well as 



the reversibility of the time evolution are no longer pos- 
sible if the effective Hamiltonian deviates from the ideal 
form ((5J. This allows us to experimentally induce local- 
ization effects by concatenating short evolution periods 
under a perturbation Hamiltonian S with evolution pe- 
riods under the ideal Hamiltonian Hq. For the present 
experiments, we choose £ = Hdd, and we label the du- 
rations of the two time periods tx; and To, as shown in 
Fig. [TJ When the duration tq = tq + of each cycle 
is short compared to the inverse of the dipolar couplings 
dij , the resulting evolution can be described by the effec- 
tive Hamiltonian 

n cS = (1 - p) Ho + p% (4) 

where the relative strength p = t^/tc of the perturba- 
tion can be controlled by adjusting the duration . Since 
the Hamiltonian Ho is generated as an effective Hamil- 
tonian, it always deviates from the ideal Hamiltonian. 
In the experiment, we compare the artificially perturbed 
evolution of H c s with the Hq evolution with its intrinsec 
errors. Note that the intrinsic errors do not produce lo- 
calization on the time scale of our experiments (see Fig. 

El). 

Taking this perturbation into account, and starting 
from thermal equilibrium, the state of the system at the 
end of N cycles is 

p n »« (Nt c ) = Up z U N , (5) 

where Un — exp | — jtH c sNtc^ is the evolution op- 
erator for the perturbed evolution. The NMR sig- 
nal, which is measured after the backward evolution 

Vjy = exp |^%o-WYo}> can ^ e wr itten as S(Nt c ) — 
Tr \Ap ndt (Nt c )\, where 

A = V N I z Vl=p n °(NT Q ) (6) 

is the effective observable and p n ° the density operator 
of the unperturbed evolution. We again determine the 
cluster size by applying rotations <p z around the z-axis, 
as in Eq. ([3]). The resulting NMR signal is then 

S(<P,Nt c )=J2^ 4 ' Ma m 

M 

= J2 e*+ M Tr {p# (Nr Q ) p%* (Nr c )} . (7) 

M 

For ideal evolution (p = 0), the individual terms Am in 
the last equation correspond to the squared amplitudes 
of density operator elements p M ° {Ntq) with coherence 
order M. For perturbed evolution, (p ^ 0), they are 
reduced by the overlap of the actual density operator el- 
ements p M dl (Ntc) with the ideal ones. To extract these 
amplitudes from the experimental data, we perform a 
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Fourier transform with respect to 4>. Two examples for 
the resulting Am are shown in the insets of Fig. [2] 

Experimental results. — Experiments were performed 
on a home-built solid state NMR spectrometer with a 1 H 
resonance frequency of 300 MHz. The spins are the pro- 
tons of polycrystalline adamantane where the strength of 
the dipolar interaction, quantified by the second moment 
of the resonance line is 7.9 kHz. In the experiments we 
chose To = 57.6/xs. The black squares of Fig. |3^ shows 
the averaged number of correlated spins as a function of 
time for an unperturbed evolution, p = 0. The observed 
cluster size K(Ntc) grows almost exponentially over the 
range considered here [27j|. The other symbols of panel 
(a) show the evolution of the number of correlated spins 
for different values of p. Initially, the cluster size K(Ntq) 
starts to grow as in the unperturbed evolution, but then 
it saturates after a time that decreases with increasing 
perturbation strength p. We consider this as evidence 
of localization due to the perturbation. The size of the 
cluster at which this saturation occurs is also determined 
by the strength of the perturbation: increasing pertur- 
bation strength reduces the limiting cluster size. Panels 
b and c of Fig. [3] visualize this localization directly by 
comparing the generation of high-order multiple quan- 
tum coherences for unperturbed (panel b) and perturbed 
(panel c; p = 0.108) evolution: they give a color-coded 
representation of the amplitudes Am (Ntc) as a function 
of evolution time Ntq- While the distribution spreads 
continuously in panel b, it reaches a limiting value in 
panel c. 

While these experiments show that the cluster size 
reaches a stationary value, they leave open the ques- 
tion if this limiting size results from a slow-down in the 
growth [21| or it represents a dynamic equilibrium state. 
We therefore repeated the above experiment for a series 
of initial conditions corresponding to different clusters 
sizes. Figure |4^i shows the corresponding pulse sequence: 
The initial state preparation, consisting of an evolution of 
duration NqTq under the unperturbed Hamiltonian Ho, 
generates clusters of size Kq. During the subsequent per- 
turbed evolution of duration Ntc, these initial clusters 
grow or shrink. Figure [4}d shows the results for two per- 
turbation strengths, p = 0.034 and p = 0.065. The filled 
symbols correspond to uncorrelated initial states and the 
empty symbols to various initial cluster sizes Kq. The 
experimental results clearly show that, for a given per- 
turbation strength, the size of the spin clusters tends 
towards the same limiting value, independent of the ini- 
tial condition. We verified this behavior for additional 
perturbation strengths (data not shown in the figure). 

Figures [3^ and0j5 indicate that the size of the resulting 
clusters decreases with increasing strength of the pertur- 
bation. To establish this dependence in a quantitative 
manner, we determined the size of the localized clusters 
from the data shown in Fig. [3] and plotted them against 
the perturbation strength (black squares in Fig. [5]). The 
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Figure 3: (Color online) a) Time evolution of the cluster size. 
The black squares represent the unperturbed time evolution 
and the other symbols correspond to different perturbation 
strengths according to the legend. b,c) Distributions of the 
amplitudes Am (Ntc) for unperturbed dynamics (b; p = 0) 
and a perturbed evolution (c; p — 0.108) respectively. The 
perturbed evolution in panel c shows localization at a cluster 
size _Kj oc ~56 spins. 
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Figure 4: (Color online) a) NMR pulse sequence for prepar- 
ing different initial clusters sizes and subsequently evolving 
them in the presence of a perturbation, b) Time evolution 
of the correlated cluster size starting from different initial 
sates. Filled symbols are evolutions from an uncorrelated ini- 
tial state for two different perturbation strengths given in the 
legend. Empty symbols start from an initial state with Kq 
correlated spins. 
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Figure 5: (Color online) Localized cluster size -Ki oc (square 
symbols) of correlated spins versus the perturbation strength 
p. Three dynamical regimes for the evolution of the clus- 
ter size are identified depending of the number of correlated 
spins compared with the perturbation dependent localization 
value: a) a cluster size decreases, b) a cluster size increases c) 
stationary regime. 



diagonal line in Figure [5] represents a linear fit to the 
experimental data represented by the black squares; its 
width indicates the error of the fit. A functional depen- 



dence K] 



loc ~ P 



-1.86±0.05 



is obtained, indicating that the 



size of the localized clusters decreases with the square of 
the perturbation strength. The limiting value for p = 1 
is K\ oc ~ 1, indicating that the system becomes com- 
pletely localized if the perturbation strength is signifi- 
cantly larger than the unperturbed Hamiltonian. The 
figure also summarizes the evolution of the cluster size 
before the static (localized) size is reached: If the initial 
size is larger than the stationary value for the given per- 
turbation strength, Kq > Ky oc , the cluster shrinks (in- 
set a in the figure, above the diagonal). If it is smaller, 
Kq < Ky oc , the size increases (inset b, below the diago- 
nal). 

Discussion and Conclusions. — Decoherence has long 
been recognized to limit the time for which quantum in- 
formation can be used. Here, we have shown that it also 
limits the distance over which quantum information can 
be transferred. To demonstrate this effect, we have com- 
pared the spreading of information in a system of nuclear 
spins under the influence of a Hamiltonian that transfers 
information and a perturbation Hamiltonian of variable 
strength. In combination, these opposing forces result in 
a quantum state that becomes localized. The localization 
size decreases with increasing strength of the perturba- 
tion. Our experimental result of a dynamic equilibrium 
size of the localized state differs from theoretical pre- 
dictions that only indicate a slow down of the spreading 
21|. The experiments were performed with nuclear spins, 
which we use as a quantum simulator, and the perturba- 
tion is taken as a model for the disorder considered in 
the discussion of localization 119-231. 



These results may also be connected to our earlier find- 
ings that the decoherence rate of quantum states with 
many correlated qubits increases with the size of the sys- 
tem [8j , indicating that larger systems are more sensitive 
to perturbations. As the system size increases, the ten- 
dency for the system to spread is therefore balanced by 
the restriction due to the perturbation. As a heuristic 
argument, we note that in a suitable interaction repre- 
sentation, the perturbation will cause a decay whose rate 
may be calculated by second order perturbation theory. 
We expect there a quadratic dependence on the pertur- 
bation strength that could be the source of the dynamic 
equilibrium size behaviour. The results presented here 
provide information about the spatial bounds for trans- 
ferring quantum information in a spin network and indi- 
cate how precise manipulations of large quantum systems 
have to be. 
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